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1 .  ABSTRACT 
1 . 1  Objective 

The  objective  of  this  effort  has  been  to  provide  to 
image  processing  systems,  such  as  the  RADC  AFES  (Automatic 
Feature  Extraction  System)  and  MIES  (Multi-Imagery 
Exploitation  System)  testbeds  an  increased  capability  for 
discriminating  tactical  targets  from  the  surface  background 
in  thermal  infrared  (IR)  imagery.  The  approach  has  been  to 
investigate  the  present  theories  of  the  thermal  emissivity 
and  reflectivity  of  natural  surfaces  and  to  apply  these 
theories  to  the  development  of  a  background  model.  This 
will  hopefully  lead  to  the  construction  of  algorithms  to  be 
used  as  a  background  information  filter  to  aid  in  the 
separation  of  targets  from  background. 


1  . 2  Background 


Target  detection  for  thermal  IR  imagery  is  usually 
accomplished  by  pixel  intensity  thresholding  routines.  net 
targets  (vehicles  with  engines  running)  appear  as  high 
intensity  (bright)  areas  on  IR  imagery.  A  threshold  limited 
detector  shows  all  pixels  above  a  specified  intensity  value 

as  targets  that  have,  at  the  time  of  imaging,  a  large 

* 

temperature  differential  from  the  background.  However,  when 
a  significant  temperature  differential  is  not  present 
targets  are  very  difficult  to  detect  by  simple  thresholding. 


ocone 


This  effort  has  addressed  the  problem  of  detecting 
targets  which  are  close  to  ambient  temperature.  This  was 
done  by  modelling  the  background  and  by  using  the  model  as  a 
base  for  developing  an  understanding  of  physical  processes 
and  their  variations,  which  will. (with  further  work)  lead  to 
the  construction  of  a  filtering  algorithm.  This  algorithm 
would  be  used  to  eliminate  the  background  information 
causing  any  objects  or  targets  with  characteristics 
different  from  the  model  to  stand  out.  We: 

a.  Investigated  theories  of  emissivity  and  reflectivity  of 
natural  surfaces  and  determined  which  would  best  fulfill 
this  effort's  objective. 

b.  Adapted  chosen  theory  to  an  imagery  format. 

c.  Developed  a  simple  model  which  simulates  the  emissivity 

and  reflectivity  of  natural  surfaces,  leading  to  an 
understanding  of  the  physical  factors  (and  their 

variations)  controlling  the  remote  sensing  process. 

d.  Developed  an  algorithm  for  the  model,  to  study  major 
remote  sensing  variables  and  their  random  and  systematic 
variations,  so  that  advances  may  be  made  towards 
developing  an  optimum  parametric  envelope  for  data 
acquisition  and  analysis. 


e.  Made  recommendations  for  further  relevant  wcr.-;. 


f.  In  this  initial 

examination 

of 

the 

problem  we 

were 

forced  by  lack 

of 

data 

and 

by 

funding  and 

time 

constraints  to 

make 

the 

following 

assumptions 

for 

numerical  calculations: 

(i)  goniometric  isotropy  of  the  target  and  background 
emissivities . 

(ii)  spectral  invariance  of  tne  emissivity  of  ground  target 
and  background  materials. 

(.iii)  transmission  through  a  homogeneous  atmosphere,  using 
typical  published  data. 

(iv)  zero  atmospheric  self-radiance 


Further  work  is  needed  in  which  improved  simulation  models 
are  used.  Indeed,  an  end-to-end  simulation  of  the  remote 
sensing  process  is  needed.  It  is  necessary  to  consider  the 
interaction  of  sensor  point  spread  function  with  tne 
heterogeneous  field  of  view:  something  which  was  not 


possible  in  this  project. 
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al  processing  can 
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for  target  detection, 
possibly  for  target 
It  is  a  technique  which 

the  target  differs  by 
from  the  background  and 
and  the  background  show 
approach  or  exceed  the 
of  the  temperatures  of 
Under  such  situations, 
)  any  other  form  of 
be  useful  in  contrast 


One  purpose  of  this  study  has  been  to  concisely  and 
clearly  state  the  established  principles  which  govern  the 
detectability  of  different  targets  from  each  other  and  from 
their  backgrounds  in  order  to  better  understand  factors 
determining  the  parametric  envelope  for  optimum  data 
acquisition  and  extraction.  The  reasons  for  this  approach 
will  become  clear  later  in  the  report. 


Those  factors  controlling  the  remote  sens 

* 

all, need  to  be  considered  in  concert,  since  the 
sensor  output  obtained  when  viewing  the  target 
to  viewing  the  background  will  depend  upon  the 
each  of  those  parameters  contributing  towards 


ing  process 
variance  in 
as  compared 
variance  in 
the  sensor 


output.  These  are,  for  example,  target  and  background 
emissivity  and  their  variances  with  angular  regime  and  with 
slope,  aspect,  soil  moisture,  microclimate  and  with  a 
variety  of  other  environmental  factors.  Also  important  are 
the  variability  of  atmospheric  transmission  in  the  bandpass 
or  combination  of  bandpasses  used  to  observe  the  target  and 
the  background;  transmission  (and  atmospheric  self  radiance) 
including  the  presence  of  scene  obscurants  .  The  interaction 
of  the  spectral  radiance  from  the  target  with  the  spectral 
response  of  the  sensor  will  be  important  in  determining 
sensor  output.  The  polarization  of  the  radiance  from  the 
target  and  its  interaction  with  the  atmosphere  (depolarizing 
effect  due  to  multiple  scattering)  and  with  the  inherent 
polarization  of  the  sensor  will  be  important  in  determining 
sensor  output.  Sky  radiance  will  be  an  additional  factor  in 
determining  composite  sensor  signal  output. 

In  the  case  of  a  heterogeneous  scene,  consisting  of  a 
number  of  scene  elements,  each  of  which  has  its  own  angular 
anisotropy  in  emissivity  and  in  bidirectional  reflectance 
(optical  reflective  regime)  the  point  spread  function  of  the 
detector  or  of  each  element  of  the  detector  (in  the  case  of 
an  array)  will  interact  with  the  heterogeneity  of  the  target 
and  with  the  angular  disposition  from  which  it  (the  target) 


is  viewed. 


Sensor  ncise,  optical  effects  and  stray  radiation 
within  the  sensing  device,  tracking  errors,  the  effects  of 
SKy  radiance,  cross-talk  between  different  elements  of  an 
array,  the  magnitude  and  dimensions  of  the  point  spread 
function  and  its  interaction  with  the  composite  scene, 
consisting  of  a  number  of  different  scene  elements  will  all 
effect  the  overall  output  signal  obtained  for  each  pixel 
viewed  by  each  instantaneous  field  of  view,  and  will 
therefore  affect  the  accuracy  with  which  target  and 
background  may  be  discriminated. 

Due  to  the  short  duration  of  this  contract,  it  was 
considered  important  to  consider  those  elements  most 
important  to  the  development  of  the  ideas  prevalent,  indeed 
predominant  in  the  statement  of  work.  To  this  end, 
classical  blackbody  radiation  theory,  published  atmospheric 
transmission  data  and  multivariate  statistical  analysis  were 
used  to  determine  the  relative  advantages  of  single  band  and 
two-band  data  acquisition  techniques  for  target  detection, 
tracking  and  quantification.  To  the  extent  that  detection 
implies  discrimination  one  two  terms  in  this  report  w'ill  be 
considered  to  be  synonymous.  The  effects  of  atmospheric 
se lf-radiance  and  the  effects  of  the  point  spread  function 
were  not  considered  for  reasons  of  time. 


Mathematical  analysis  has  been  performed  for  single 
band  and  for  band  ratio  techniques  in  order  to  determine  the 
feasibility  of  target  discrimination  at  low  target-to- 
background  temperature  differences.  The  examples  selected 
are  specific  to  a  limited  range  of  target  and  background 
emissivities ,  which  are  considered  to  be  Lambertian. 
Atmospheric  transmission  is  considered  to  be  uniform  and  to 
consist  of  certain  pre-selected  values  taken  from  the 
literature  as  typical.  Variances  for  these  parameters  have 
been  selected  from  an  examination  of  the  literature  (e.g. 
the  Infrared  Handbook  and  are  considered  to  be  reasonable 
for  a  first  order  analysis. 

A  literature  survey  of  the  current  state  oi  the  art  as 
reported  in  the  unclassified  literature  has  been  performed 

and  is  included  in  the  discussion.  Conclusions  may  be  drawn 
with  respect  to  promising  new  avenues  of  research  for 
improved  target  tracking  and  discrimination,  for  the 
development  of  an  optimum  parametric  envelope  for  data 
acquisition  and  analysis,  for  optimized  target  detection 
tracking  and  quantification.  Suggestions  are  made  for 
continuation  and  extension  of  this  work  in  related  areas 
which  we  had  jao  time  to  investigate  in  this  study  such  as, 
for  example,  dynamic  viewing  using  staring  two-dimensional 
pyro- electric  or  quantum  effect  arrays  using  CID  or  CCD 
devices.  While  the  reported  literature  shows  empirical 
comparison  and  evaluation  of  dynamic  viewing  of  selected 
targets  under  ideal  conditions,  it  is  considered  that  far 


more  extensive  modelling  studies  are  required  in  oru-_-r  re 
determine  the  optimum  parametric  envelope  for  using  two- 
dimensional  arrays  in  a  dynamic  viewing  method  of  operation. 

These  remarks  apply  to  nadir-pointing  sensors,  but  will 
apply  even  more  to  reasonably  wide  angle  (10  -  20  degree 

field  of  view)  scanned  or  two-dimensional  sensor  systems 
viewing  ot-_quely  towards  (say)  the  horizon  from  space  or 
from  airborne  platforms  ,  or  even  at  ground  level. 

A  major  emphasis  in  this  report  is  to  point  the  way 
towards  future  areas  in  which  it  is  recommended  that 
research  efforts  be  concentrated. 

A .  FACTORS  CONTROLLING  THE  RELATIONSHIP  OF  REMOTELY  SENSED 
RADIANCE  TO  GROUND  FEATURES 

A . 0  Introduction 

While  there  has  been  a  consideration  of  the  many 

factors  affecting  the  remote  sensing  process  (e.g.  ^  )  and  a 

consideration  of  those  factors  affecting  the  selection  of 
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:  impasses  for  smart  sensors  *  ’  there  has  been  little 

4-9  ,  -  10,1 

wcrK  on.  tne  effects  of  systematic  and  or  random 

variations  on  the  accuracy  with  which  targets  may  be 
discriminated  using  radiance  data.  We  sha 1 1 - cons ide r  here 
each  of  the  factors  controlling  remotely  sensed  radiance 
levels,  the  equations  describing  those  factors  and  the  level 
of  variation  which  may  be  typically  anticipated  to  occur  in 


c  -J 

tally  e 

f  i  e  c  c  1  v  . 

huge  r. ,  Ten 

o c  n  ana  Gray 

snowea  t.nat 

up 

to  30 % 

of  the 

IF  0V 

cou  la 

be  filled  by 

cloud  over  a 

-y 

pica  j. 

target, 

before 

the 

resulting  pixe-  o-: 

recognized  as  cloud,  using  current  NOAA  AVHRR  da  tc  screening 
algorithms.  Cloud  contaminated  pixels,  once  recognized,  are 
discarded  from  further  analysis,  since  there  is  no  accurate 
way  to  correct  for  cloud  without  making  assumptions  about 
tne  degree  of  cloud  contamination  of  flagged  pixels.  A 
qualitative  illustration  of  the  effects  of  cloud 

contamination  is  to  point  out  that  areas  adjacent  to  cloud- 
covered  regions  usually  have  a  modified  color  balance  on 
interactive  computer  screens.  That  is,  such  areas  may  be 
viewed  as  "different".  However,  a  far  more  quantitative 
test  than  human,  viewing  of  an  interactive  computer  screen  is 
u-.-d  to  ueiete  pixels  which  are  cloud  contaminated,  or  are 
suspected  of  being  so. 

investigations  are  currently  underway  to  study  the  of  f- 
nadir  effects  of  unresolved  cloud  and  haze  on  target 
discrimination  for  atmospheres  of  different  turbidity. 


.  1  .  .  Tne  effect  of  random  vs.  systematic  errors  in 

r_e  corded  radiance  levels. 

6  6 

Ac  pointed  out  by  Duggin  and  by  Duggin,  et.  al 

t h-.-re  arc  both  systematic  and  random  errors  in  radiance 
levels  recorded  from  targets  which  are  to  be  discriminated 

and  cuan  t : :  i  <•  d.  The  systematic  effects  may  be  corrected  for 


This  follows  ircx  the  general  case  of  a  mixed 
(heterogeneous)  pixel,  as  shown  in  Table  1.  Here  is  the 
traction  of  the  instantaneous  field  of  view  (IFOV)  filled  by 
target  A,  while  a  is  the  fraction  of  the  IFOV  filled  by 
cloud.  Thus,  in  addition  to  atmospheric  effects  mentioned 
previously,  opaque  cloud  (or  haze  or  cirrus)  can  alter  the 
absolute  radiance  ievei  recorded  when  viewing  a  given  target 
and  may  (depending  on  tne  spectral  reflectance 
characteristics  of  the  cloud)  alter  the  spectral 
distribution  of  target  radiance. 

Cloud  which  fills  only  a  portion  of  the  IFOV  of  the 
detector  does  not  result  in  a  radiance  level  or  in  a 
spectral  distribution  in  radiance  recorded  by  several 
channels,  which  characterizes  cloud.  However,  unresolved 
cloud  can  have  the  same  effect  as  cirrus  or  haze,  namely  to 
distort  the  spectral  radiance  signature  from  the  target, 
which  is  to  be  discriminated  from  its  surroundings  or 
quantified  in  some  manner.  While  the  photointerpreter 
relies  upon  shapes,  textures,  relative  positions  and 
i.  oo  tii  or. s  as  we  1 1  as  on  color  and  brightness  ,  computer 
class:  *'i;eticr  of  mul  tispectral  imagery  relies  only  upon  the 
digital  counts  recorded  from  the  target  in  each  sensor 
channel.  Therefore,  unresolved  cloud,  cirrus  and  haze  can 
result  in  errors  of  emission  and/or  of  commission  if  there 
is  not  some  means  of  detecting  the  cloud,  haze  or  cirrus- 
contaminated  pixels  so  that  they  are  not  considered  in  the 
analysis.  Co  far,  cloud  s  iroening  algorithms  are  not 


oa.T.e  areas  ina  will  permit  a  amp  a  r  ..sons  to  be  sios  Letweer. 
images  obtained  over  different  areas.  Such  corrections  are 
necessary  to  permit  analysis  for  monitoring  change  in  ground 
cover  over  large  geographic  areas  and  will  enable 
quantitative  evaluation  of  temporal  change  in  the  values  and 
in  the  spectral  distribution  of  recorded  radiance. 

Such  angle-dependent  effects  are  systematic  and  may,  to 
some  extent,  be  corrected  for.  This  may  be  achieved  by 
parallel  empirical  and  tneoretical  (modelling)  studies. 


B.I.e.  The  effect  of  cloud  and  haze. 


Cloud,  haze  and  cirrus  may  be  present  in  constant  or  in 
predictably  varying  quantities  across  a  scene,  in  which  case 
they  will  modify  recorded  radiance  as  reported  by  Duggin, 
et.  al  ^  ^  .  The  overall  recorded  radiance  signature  in 
bandpass  r  will  be,  from  a  given  target  (A)  (given  that  the 
cloud  is  a  Lambertian  reflector  with  a  spectrally  invar¬ 
iant  reflectance  factor  R^) . 
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(3.11).  1(a)  is  the  spectral  response  of  the  sensing  device 
and  equation  (3.12)  may  be  more  fully  expressed  (so  that  the 
relative  importance  of  the  abovementioned  factors  may  be 
appreciated)  as: 


1(A). dA 


(3.13) 


Here  the  sensor  view  and  azimuth  angles  are  6'  and  e '  . 

Not  surprisingly,  it  has  been  shown  that  there  is  a 

significant  dependence  of  recorded  digital  radiance  levels 

on  sun-target-sensor  geometry  for  the  NOAA  AVHRR,  SPOT 

A  63 

and  MSS  data  .  There  are  indications  that  this  is  also 
true  for  thematic  mapper  (TM)  data.  Empirical  understanding 
of  the  dependence  of  recorded  radiance  on  scan  angle,  season 
(solar  declination)  and  latitude  should  lead  to  an 
approximate  correction  for  these  parameters,  permitting 
radiometric  comparisons  (and  therefore  classifications  on 
the  basis  of  radiance  levels)  across  images.  Such 
corrections  will  also  permit  quantitative  comparisons  to  be 
made  between  images  obtained  at  different  seasons  over  the 


There  are  "zero  power  wavelengths"  (1,  ,  .<2  )  which  bound  each 
sensor  bandpass,  and  at  which  the  response  of  the  sensor  is 


zero.  However,  between  these  wavelength  limits,  the 
response  of  the  detector  is  not  constant.  Since  the  target 
radiance  and  the  atmospheric  transmission  and  path  radiance 
also  vary  (each  in  a  different  manner)  across  each  bandpass, 
interactions  between  the  spectral  response  of  the  sensor, 
the  spectral  radiance  from  the  target,  the  spectral 
backscatter  from  the  atmosphere  and  the  spectral  extinction  of 
the  atmosphere  along  the  path  from  the  target  to  the  sensor 
may  be  anticipated.  This  has  been  pointed  out  by  Slater^ 
and  by  Duggin,  et.  al  10,ol,6z.. 


The  signal  output  ny  the  sensor  (normalized  sensor 
response)  is  given  by  the  equation: 


(3.12) 


where  L  (  1  ,  the  spectral  radiance  incident  on  the 
sensor  is  due  to  radiance  reflected  from  the  target  for  the 
sun- target-sensor  geometry  considered,  after  modification 
for  the  spectral  extinction  along  the  atmospheric  path  from 
the  target  to  the  sensor,  plus  the  added  spectral 
atmospheric  backscatter  (path  radiance)  for  the  same  target- 
sensor-geometry  as  considered  in  equations  (3.10)  and 
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wnere  -  '  v  ",  :  ",*, z)  is  the  extinction  coefficient  j- 1 

fci  X.  L 

altitude  Z  along  the  path  to  the  sensor  and  L  ; 6  1  ,$■  1  , A ) 

is  the  path  radiance  scattered  into  the  sensor  by  the 
atmosphere  in  addition  to  radiance  transmitted  from  tr.e 
target  at  ground  level. 

It  may  be  readily  seen  that  the  contrast  in  recoroec 
radiance  between  two  targets  is  reduced  by  botn  a tmospnc r ± c 
extinction  and  by  the  addition  of  a  background  radiance 
level  (e.g.  ref.50  ,  pg.  206).  The  reason  for  this  is  that 
the  level  of  useful  signal  is  reduced  by  the  atmosphere, 
while  there  is  an  addition  of  purely  atmospheric  signal 
which  contains  no  useful  information. 

Some  correction  can  be  made  for  path  radiance  and  ior 
atmospheric  extinction.  However,  methods  currently  use_ 
assume  that  the  atmosphere  is  constant  across  the  imaged 
area  (e.g.  ^5).  Further,  either  approximations  are  made  in 
that  all  signal  from  low  albedo  areas  is  supposed  to  arise 
from  the  atmosphere  (e.g.  ^5  )  or  several  meteorological 

inputs  are  required  in  order  to  compute  the  modelled 
corrections  (  e  .  g  .  5  )  . 
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B . 1 . d .  The  effect  of  the  detector. 

The  response  of  the  detector  is  wavelength- 
dependent  across  those  bandpasses  to  which  it  is  sensitive. 
In  other  words,  it  is  not  either  "on"  or  "off",  depending 
upon  wavelength.  This  is  illustrated  in  Figs.  9  and  1C. 
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gor.ioir.-i trie  nemispher ical-directional  reflectance  factors  ' 
Variations  in  global  irradiance  across  an  imaged  area 
will  also  occur  due  atmospheric  transmission  fluctuations 
(tnat  is,  variations  in  time  at  a  given  location  will  be 
similar  to  variations  with  space  at  a  given  instant  in  time) 


Radiance  froc 


the  ground  will  also  be  modified  by  the 


atmosphere,  prior  to  detection. 


B.I.c.  The  passage  of  reflected  target  radiance  to  the 


detector 


Radiance  reflected  from  the  target  may  be  most  exactly 

expressed  by  equation  (3.10).  However,  some  radiance  is 

absorbed  and  scattered  out  of  tne  reflected  beam  by  the 
77  -  128 

atmosphere  .  This  rnay  be  considered  by  using  a 

modification  of  equation  (3-2)  allowing  for  the  fact  that 
extinction  along  the  reflected  path  may  differ  from 
extinction  along  the  path  of  the  incident  radiance  falling 
on  the  target.  Fig.  8  explains  this  relationship.  Thus,  the 
radiance  falling  on  the  detector  is  given  by: 
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The  relationship  between  incident,  reflected  and  absorbed 


radiance  is  shown  in  Fig.  7.  As  pointed  out  by,  for  example, 
19  52 

Kriebel  and  by  Kirchner,  et  al  ,  it  is  the  interaction 

of  the  polar  distribution  of  the  spectral  irradiance  field 
with  the  bidirectional  spectral  reflectance  factor  which  deter 
mines  the  angular  distribution  (and  the  absolute  value  at  any 
given  angular  orientation)  of  spectral  reflected  radiance. 
Thus,  equation  (3.9)  may  be  more  exactly  expressed  as: 
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where  R ( 0  ,  4> ;  0 


+  LT(e,*,X).R(Q,*;e',*',x) 

c  o 

<$>')  is  the  bidirectional  reflectance  factor J 


Thus,  while  equation  (3.9)  describes  the  spectral 
radiance  reflected  into  a  given  direction,  as  a  function  of 
overall  irradiance  level,  equation  (3-10)  takes  into 
consideration  the  polar  distribution  of  the  incident 
spectral  irradiant  flux  field. 

It  has  been  shown  (e.g.  10*11»19,52  ^  that  the  scalar 

global  spectral  irradiance  and  its  polar  distribution  are 

time-dependent.  Therefore,  both  the  scalar  value  and  the 

* 

polar  distribution  of  this  quantity  may  vary  between 
sequential  ground  measurements  of  target  radiance  and 
irradiance,  giving  errors  in  estimates  of  spectral 


Tne  radiance  reflected  into  any  polar  direction  may  be 
expressed  as: 


L  (9',<>\A)=R(2ir  E(6,*,X) 


where  R( 2  H;  0'  ,  b'  X)  is  the  hemispherical-directional  spectral 
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reflectance  factor  ,  which  relates  the  spectral  radiance 
reflected  into  any  selected  direction  (selected  oy  tne  scan 

angle  and  by  the  ground  track  of  the  sensor  platform)  to  the 
total  scalar  (global)  spectral  irradiance.  The  angular  noca 
tion  used  is  shown  in  Figs.  A  and  5. 


There  can  be  considerable  angular  anisotropy  in  the 
spectral  hemispherical-directional  reflectance  factor. 
Indeed,  a  polar  plot  of  this  factor  generally  produces  a 
non-spherica 1  surface,  whose  irregularity  is  indicative  of 
the  degree  of  angular  anisotropy  in  the  spectral 
bidirectional  reflectance  factor.  Experimental  evidence  of 
such  anisotropy  has  already  been  widely  published  (e.g. 
furthermore  ,  theoretical  calculations  indeed  predict  such 
anisotropy  (e.g.  )  in  a  manner  which  has  been 
experimentally  confirmed  (e.g.  ^  ).  Fig.  6  shows  the 
variation  of  the  bidirectional  reflectance  factor  with 
wavelength  and  scan  angle  for  wheat  at  a  growth  stage  of  3. 5 
on  the  modified  Feeks  scale,  calculated  from  published  data 
^  .  The  solar  azimuth  and  zenith  angles  are  taken  to 
correspond  to  those  for  the  NOAA-6  advanced  very  high 
resolution  radiometer  (AVHRR). 


wiiere 


h  --  a ’i  i  i  cude  (km) 

2 

<s  r  *  Knyleigh  scattering  cross  section  (ra  ) 

-  3 

nr  -  acmospneric  number  density  (m  ) 
np  *  aerosol  number  density  (m~^) 

Ay  =  ozone  absorption  coefficient  (cm-1) 

D-j  =  ozone  concentration  (cm  km"^)  . 


B.I.b.  The  reflection  of  radiance  from  the  target. 

The  situation  at  the  target  is  shown  in  Fig.  3-  The 
total  resolved  component  of  radiance  falling  perpendicularly 
on  the  surface  is  shown  as  the  scalar  quantity  E  ( X) .  This 
quantity  will  depend  upon  the  solar  zenith  angle  ,  upon  the 
various  atmospheric  extinction  coefficients  and  upon  the 
amount  and  polar  distribution  of  cloud  (  ■  )• 

There  will  be  some  energy  which  will  be  absorbed  by  the 
ground  target.  This  is  E  ^  (  8 '  , 4> '  ,  \  )  and  will  depend  upon 
polar  direction,  since  there  will  probably  be  angular 
anisotropy  in  the  absorbing  properties  of  the  ground. 

The  energy  reflected  in  a  given  polar  direction  will  be 
L  ^  (  e'  ,  4>'  ,  \)r  There  will  be  a  net  balance  as  given  by  the 
equation: 

n  n  n  « 
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W1'ere  Sext(6"’:"’A,2)  *  sr(eM.«";x,*)  +  6p(6".<r ,x.i)  +  63(e".0''^,2) 

Here  (  0"  .  a  •  2)  is  Che  total  extinction  coefficient 

(8"»*£".A,z)  is  the  coefficient  for  Rayleigh  scattering. 

6p  A,z)  is  the  coefficient  for  aerosol  scattering. 

(0".$",X,2)  is  the  coefficient  for  ozone  absorption. 

The  total  irrradiance  on  the  target  at  wavelength  X  is  given 
by  the  sum  of  the  radiance  from  the  sun's  disk  L^(Q  ,<p  ,X) 
and  the  irradiance  from  the  sky  (scattered  radiation),  see 


(3.3) 


Figs.  2  and  3. 


K  (X)  -  Lt  (0,<M>  Cos  0  + 


ffl 

/  /  scatt 


(e,4i,A)  Cos  0  d  0.d  <p 


(3.4) 
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Here  E  (  X)  is  a  scalar  quantity;  while  E  scatt  (6,$,X)  is  a 
polar  vector  which  represents  the  scattered  radiance  from 
each  point  in  the  sky.  Integration  of  the  component  Escatt^  ’  ^ 
over  the  hemisphere  of  the  sky  gives  a  scalar  quantity 
which  represents  the  total  "sky"  irradiance.  L^,(  6,  <f>,  A)  .Cos  6 
is  the  resolved  component  of  the  radiance  from  the  solar 
disk  which  falls  perpendicular  to  the  ground. 

Approximate  values  for  the  extinction  coefficients,  in 
terms  of  readily  measured  physical  quantities  are  given  by 
(e.g.)  Slater_"^. 


3  r  (h)  =*  a  rnr(h)  x  10' 
a  p(h)  =  ^  p(°> 


(3.5) 


(3.6) 


0  ■»  ( h  )  ~  A  D  j  (  h ) 
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5.1  The  optical-reflective  region  of  the  spectrum 


In  the  optical-reflective  region,  we  shall  consider  the 
remote  sensing  process  in  three  parts;  the  passage  of 
radiant  energy  from  the  sun  to  the  target,  the  reflection  of 

energy  by  the  target  and  the  interaction  of  the  atmosphere 
with  reflected  radiance  before  it  reaches  the  sensor. 


B.I.a.  The  passage  of  radiant  energy  to  the  target. 


Fig.  2  shows  the  passage  of  solar  radiance  through  the 
atmosphere  to  the  target  on  the  ground.  The  solar  zenith 
and  azimuth  angles  are  determined  by  the  local  time,  season 
(solar  declination  d)  and  latitude  C  .  In  fact,  the  solar 
zenith  angle  6  is  given  by  the  equation: 

Cos  6  =  Sin  *.Sin  6  +  Cos  <J>.Cos  6  Cos  15  (TH  +  ^  -  SLN)  (3.1) 


where  TH  is  the  local  time  in  hours,  TM  the  local  time  in  minutes  and 
where  SLN  is  solar  noon,  the  time  at  which  the  sun  crosses 
the  meridian. 


Losses  from  the  radiant  beam  from  the  sun's  disk  occur 
due  to  scattering  and  absorption.  Broadly,  the  radiance 
reaching  the  target  along  the  direct  beam  from  the  sun's 
disk  is  (e.g.  ^). 


L.r(e , : ,  a) 


Lo(0.*.*) 


,  x.e)  .de".d<r.  dz 


(3.2) 


eilecc  of  the  atmosphere,  in  order  tc  correct  for  it  ( e .  g . 
i/’  /-'4  ).  Such  models  generally  assume  that  the  atmosphere 
is  not  turbulent,  but  rather  consists  of  a  series  of  steady 
state,  layered  media,  with  discrete  and  constant  aerosol, 
molecular  and  particle  concentrations. 


Before  dealing  with  variations  in  these  quantities,  we 
must  quantitatively  describe  the  entire  remote  sensing 
process  in  an  overall  manner  as  outlined  by  (e.g.)  Maxwell1. 


ome  of  those  factors  (on  the  basic  of  what  very  limited 
experimental  data  exists).  The  effect  of  such  variation  on 
recorded  radiance  levels,  on  target  discrimination  and 
quantification,  on  system  design  and  on  data  collection  and 
analysis  optimization  will  be  considered. 

B . 0  Discussion 

The  interactions  occurring  between  solar  radiance 

incident  on  the  top  of  the  atmosphere,  the  atmosphere,  the 

ground  and  the  sensor  are  shown  broadly  in  figure  1.  It  is 

seen  that  solar  radiance  ( L  )  incident  on  the  top  of  the 

earth's  atmosphere  is  subject  to  losses  in  each  polar 

direction  (  5  ,  4>  )  caused  by  scattering  (Lg  (6  ,  :  )  )  ana 

13-19 

absorption  (L  (6,4:))  (e.g.  )  during  its  passage  through 

the  atmosphere.  The  total  energy  which  reaches  the  earth's 
surface  is  given  by  the  sum  of  the  radiance  from  the  solar 
disk,  plus  the  energy  scattered  from  the  sky  to  the  earth's 
surface  (diffuse  irradiance  or  "sky  radiance").  The  earth's 
surface  reflects  radiation  in  an  anisotropic  manner,  which 
has  been  described  analytically  20-25  ancj  which  has  been 

0/1  _  o/l 

observed  experimentally  (e.g.  Q  ”  J  ).  The  radiance 

reflected  from  a  ground  target  into  the  direction  of  the 
sensor  is  again  subject  to  scattering  and  absorption.  Some 

<e 

radiance  reflected  from  adjacent  target  areas  is  scattered 
into  the  reflected  beam  and  radiance  is  scattered  from  the 
atmosphere  into  the  sensor  ("path  radiance"  or 
"backscatter" ) .  Models  have  been  developed  to  describe  the 


etnp  .riuaily ,  using  regression  equations  developed  o.. 
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studying  a  large  volume  of  data  ’  .  However,  the  nature  of 

random  errors  prevents  their  correction.  The  best  that  may 
be  achieved  is  to  determine  the  magnitude  of  trie  random 
errors  and  to  dstermine  their  effect  upon  feature 
identification  and  quantification  Tne  effects  on  sensor 

output  of  random  and  systematic  errors  in  various  remote 
sensing  variables  are  shown  in  flow-chart  form  for  tne 
optical-reflective  part  of  the  spectrum  in  Fig.  11  and  for 
the  thermal  infrared  region  in  Fig.  12. 


Thus,  in  the  event  that  two  targets,  A  and  E  are  to  be 
discriminated  and  there  is  no  visible  cloud  in  the  sky,  then 


the  radiance  difference  needed  to  perform  a  discrimination 
with  9595  confidence  would  be  given  by  the  expression 


where  H.W.  denotes  half  the  width  of  the  95$  confidence 

interval  on  the  difference  (L  )  -  (L  )  and  where  { Lr )  and 

r  A  r  B  A 

(L  )  are  the  mean  radiance  values  in  bandpass  r  fr:r. 
r  3 

targets  A  and  B.  The  means  are  obtained  in  each  case  from  n 
pixels.  ^  o.  5  m_q  is  the  student  factor  where  (m)  is  the 
number  of  observations  used  in  estimating  (  L^. )  ^  and  ( Lj.  )g  . 
The  evaluation  of  equation  (5-15)  using  equation  (5- 14) 
would  be  complex.  However,  if  we  use  the  concept  of 
hemispherical-directional  reflectance,  in  equation  (5-9) 
then  we  can  obtain  a  simple  solution  from  the  works  of 


,<J)  )  describes  the  sun-target-sensor  geometry  for  the  group  of  pixels  in  target 


Here  we  refer  to  mean 


specific 
is  the 


radiance  differences  at  a 
wavelength  A,  instead  of  across  bandpass  r.  t 
atmospheric  transmission.  Note  that  we  assume  a  fraction  of 
cloud  ac  in  the  pixels  over  target  A,  while  we  assume  a 
fraction  of  cloud  a ^  in  the  pixels  over  target  E.  This 
formula  would  need  to  be  computed  for  each  scan  angle, 
putting  in  the  values  for  each  variable  relating  tc  the 
appropriate  sun-target-sensor  geometry. 

It  should  be  noted  that  if  we  consider  the  radiance 
output  from  the  sensor  (recorded  radiance  rather  than  the 
radiance  incident  on  the  sensor  as  considered  above),  then 
for  each  wavelength  ,\  equation  (3.16)  becomes,  more  fully 
expressed  in  equation  (3.19)  (next  page).  These  equations 
enable  us  to  compute  the  necessary  mean  radiance  difference 
between  two  targets  for  target  discrimination  with  a  95Z 
probability  of  being  correct.  The  computation  is  unique  to 
each  sun- target-sensor  geometry  and  will  depend  not  only  on 
the  angular  anisotropy  of  each  target  reflectance,  but  will 
else  depend  upon  the  angular  anisotropy  of  the  atmospheric 
scattering  functions.  The  computation  is  also  sensor-specific 

In  detail,  the  sensor  output  difference  A  for  m  pixels 
obtained  over  each  target,  each  partly  filled  by  a  proportion 
of  (Lambertian)  cloud  would  then  be: 


I(\).L  U)-d\  -  I  1 (X )  •  L  ( A )  •  < 

A  JXX  B 

/  l  ( X )  •  dA 

>1 


I(X)-|  La(X)  -  Lb(A)  l  -dX 


I(X)*d\ 


In  order  to  be  95%  confident  that  the  two  features  are 
indeed  seperable  on  the  basis  of  reflected  radiance  using  a 
bandpass  r,  whose  upper  and  lower  zero-power  wavelength 
limits  are  A9  and  ,  the  difference  must  exceed  the  value 
of  the  following  function: 


(<5^  ■  <r'\) 


/h.w.{ 


I(X)  .1.  (A)  -  1(A) .Lr(X)  \  .dx 


J  I(X)-< 


where 


K.W.  [(LpA  -  <Xr')B}  i 


is  the  minimum  recorded  radiance 


difference  a  necessary  to  be  95%  confident  of  discrimination 
on  the  basis  of  sensor  output.  Of  course,  we  can  perform 
such  calculations  for  any  pre-determined  confidence  level. 


Factors:  which  will  control  the  discriminabil  i  ty  of 

targets,  where 

H.w.  /(T)  -  (Tj  \ 

Discriminability  =  D.F.  =  V.  A  B  J 

function  |  I 

M 

as  given  by  (e.g.  Duggin  ^  )  are  several.  T  nese  effects 
are  mentioned  but  have  not  been  investigated  in  this  study  . 
For  discrimination  to  be  possible  D.F.  <  1.0.  ? 


(i)  relative  areas  in  the  IFOV  occupied  by  unresolved  cloud 

and  by  target. 

(ii)  atmospheric  extinction  and  its  variation  (here  shown  as 
atmospheric  transmission  t  )  for  the  sensor  bandpasses. 


(iii)  variation  in  the  relative  proportions  of  the  IFOV 
occupied  by  the  target  and  by  cloud.  ^  ^7 


' iv)  path  radiance  for  tne 
variatior.  across  tne  1:, 


•...cor  bandpasses  anc  its 
64 


(  e 


) 


(v)  relative  albedos  of  tne 


r go  -s 


fc  be  distinguished  in 


the  bahdpasses  used  and  of  the  cloud  present  in  the 
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IFOV  where  unresolved  cloud  exists. 


(vi)  variations  in  atmospheric  turbidity  and  the  average 

,  68. 

imaged  area.  \e-g*  / 


value  of  this  quantity  across  the 


.  ii.)  Tne  Overall  level  and  spectral  distribution  oz  global 
irradiance . 

(viii)  The  spectral  responses  of  tne  sensors  in  those 
bandpasses  used,  and  the  interaction  of  such  responses 

with  the  target  radiance. 

(ix)  random  variability  of  the  spectral  b idi rec t iona 1 
reflectance  factor  across  the  sampled  area.  This  can 
be  caused  by  (e.g.)  topographic  variations,  textural 
variations,  slight  changes  in  pixel  composition,  etc. 

(x)  the  dependence  on  sun-ta rge t-sensor  geometry  of  the 

spectral  bidirectional  reflectance  factor,  spectral 
path  radiance,  spectral  atmospheric  extinction  and 
target  aiscriminability .  Due  to  angular  anisotropy  of 
ground  reflectance,  it  may  be  that  discrimination  of 
some  targets  is  optimal  (target/background  contrast  is 
maximum)  at  certain  angular  regimes  (i.e.  at  certain 
solar  elevations  and  azimuths  and  at  certain,  possibly 
off-nadir  view  angles). 

(xi)  The  interaction  of  the  point-spreaa  function  of  the 
detector  with  the  heterogene  1  by  cz  zhe  pixel, 

consisting  of  different  scene  components.  This  point 
is  made  clear  in  Fig.  13*  Here,  the  point-spread 
function  of  the  sensor  is  shown  superimposed  upon  a 
nominal  pixel  (projection  of  the  defined  earlier  (IFOV) 
on  the  ground)  containing  several  different  scene 


elements .  Each  scene  element  will  nave  ice  own, 
different  angular  spectral  reflectance  and  emissivity 
anisotropy  and  so  the  movement  of  the  sensor  such  test 
the  IFOV  is  shifted  by  a  fraction  of  a  pixel  will 
produce  a  considerable  change  in  sensor  output.  This 
important  point  has  so  far  received  very  little 
attention. 

An  example  ^7  0f  the  effects  of  scan  angle, 

unresolved  cloud  and  atmospheric  turbidity  on  contrast  ratio 
is  shown  in  Figs.  14  (clear  atmosphere:  meteorological  range  > 
50  4m)  and  15  (turbid  atmosphere:  meteorological  range  <  10 

km).  Here  vegetated  targets  (70:2  wheat,  30%  soil)  and  100S 
wheat  at  a  growth  stage  of  5*5  on  the  modified  Feeks  scale 
(boot  stage)  are  considered.  In  each  case,  a  vegetative 
index  (VIN)  is  used  to  typify  the  target:  the  VIN  used  is 
the  radiance  re corded  in  AVHRR  band  2  (0.713-0.986  pm; 

reflected  infrared)  divided  by  the  radiance  recorded  in 
AVHRR  band  1  (0.570-0.686  urn;  visible).  Three  facts  emerge: 


firstly  the  contrast 

r3tic 

(VIN  /VIN 

V  MAX  'MIN 

)  is  greatly 

reduced  by 

unresolved 

cloud 

(and  therefore 

by  haze  and 

c  i  *.  1  tio  /  • 

Second  I y , 

the 

contrast  ratio  is 

scan  angle 

dependent . 

Thirdly, 

the 

contrast  ratio  is 

reduced  in 

a  scan  angle- dependent  manner  by  unresolved  cloud  for  a  turbid 
(meteorological  range  <  10  km)  atmosphere  as  compared  to  a 
clear  (meteorological  range  >  50  km)  atmosphere. 


5.2.  TH THERMAL  INFRARED  REGION  c:  THE  SPECTRUM 

Evidence  exists  to  snow  that  there  is  angular 
anisotropy  of  the  emissivity  of  ground  targets  71-76  >  in 
the  case  of  emitted  thermal  radiance,  'atmospheric 
attenuation  consists  principally  of  absorption,  so  that 
equation  (3*10)  should  be  re-written,  for  emitted  thermal 
infrared  radiance  at  ground  level,  as: 

Lr(  0  1  ,♦•  ,  A )  =  e  (0'  ,  4> 1  ,  X  )  .  W  (T ,  X  )  (3*2 

where  e  (  0  1  ,  4>  '  ,  X  )  is  the  directional  radiant  spectral 
emittance  and  W(T,  X)  is  the  spectral  radiant  exitance  of  a 

blackoody  at  temperature  T  (degrees  Kelvin).  The 

assumption  here  is  that  the  radiator  is  a  graybody. 


Therefore,  in  the  case  of  emitted  thermal  radiance,  the 
equivalent  equation  to  the  optical-reflective  equation  (3.11) 
Is  : 


L^Q’.d’.X)  -  LpCe* ,X).expl- 


6  '  ext^  3"  >  <1)"  >  ^  >  z  >  de".dd".dz 


+  L 


sky 


(0' ,X) 


(3.21) 


Equations  (3.15)  and  (3-16)  would  also  be  different,  of 
course,  for  tne  tnermal  infrared  case  than  for  the  optical 
reflective  case.  Here,  we  ignore  the  (generally  trivial) 
covariance  terms. 


HW  |LU)a  -  L ( A ) B 

■  (('*  STT  +  •»<**.*> 


£B£B 


+  (£B2  (S’.f'.A)-  STT  +  T2  (o'  ,  b '  ,  A )  •  S£  e  )  •  W(Tg , 

(3.23) 


+  2  S 


LskyLsky 


% 


where  S  ,  S  _  ,  Sr  c  and  ST  T  represent  sample 

TT  £A£A  £BeB  Lskysky 

variances  in  atmospheric  transmission,  emissivity  of  target 

A,  emissivity  of  background  B  and  sky  radiance  Lsky 

Atmospheric  transmission  t  is  given  by 

T  =  exp  f-jZ  j  J  8' ext  (9",b",A)  .d0".de".dz 
^  0  -71  0 

More  discussion  will  be  made  of  this  derivation  later. 


Equation  (3-23)  could  then  be  substituted  into 


equations  (3:18)  and  (3.19)  to  determine  the  value  of  a 
discrirainability  function  (DF)  such  that: 


35 


and  D.F.  must  always  :  .  ess  than  1  ,  for  discrimination  c: 

target  from  background  with  a  95%  (or  some  ether  pre¬ 
determined  level)  of  confidence. 


The  well-known  physical  factors  controlling  the  remote 
sensing  process  have  been  discussed  and  represented  in 
relevant  form  in  this  section.  The  relative  importance  of 
some  already  reported  random  ana  systematic  variations  on 
recorded  radiance  and  on  target  di scr iminabi 1 i ty  and 
contrast  have  also  been  discussed.  These  factors  need  to  be 
borne  in  mind  when  designing  instrumentation  or  when 
optimizing  imaging  conditions. 


Assuming  that  the  wavelengths  for  the  bandpasses  used 
cy  the  sensing  device  have  been  best  chosen  for  the  task  at 
nand,  (and  this  has  been  the  case  only  so  far  as  is  possible 
with  limited  data)  it  is  possible  that  discrimination  of  the 


targ‘ 


will  be  more  effective  for  some  sun-target-sensor 


geometries  than  for  others.  It  is  also  likely  that  there 
are  upper  limits  for  unresolved  cloud,  haze,  obscurants 
and  cirrus  for  which  discrimination  of  targets  of  different, 
selected  albedos  is  possible.  In  order  to  understand  the 
principles  controlling  both  tne  systematic  and  random 
variations  discussed  above,  which  control  the  accuracy  with 
which  selected  targets  may  be  discriminated,  it-  is  necessary 

to  use  both  empirical  and  mathematical  simulation  studies 

37-49  90-128  .  93 

for  the  atmosphere,  ’  for  sensing  devices,  ’ 

5-10.14,59-62,136,164,241-266  ,  r  .  r 

’  ’  ’  and  lor  tne  variability  of 

.  .  .  .  .  .  .  -  .  . .10-12,201-240 
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rc- 


. order  to  rest  and  calibrate  the  mathematical  models  it 
necessary  to  have  sufficient  measurements  of  the  reflectance 
factors  and  emissivity  of  the  earth's  surface  at  various  an¬ 
gular  geometries.  unfortunately,  few  such  data  currently  ex¬ 
ist  and  there  is  little  agreement  on  the  methods  by  which  such 
critical  data  should  be  obtained. 

While  it  is  necessary  to  be  pragmatic  regarding  t:.e 
best  use  of  existing  data,  it  is  suggested  that  few  advances 
in  sensor  design,  in  data  collection  optimization  or  in  an 
understanding  of  tae  effects  of  random  and  systematic  errors 
on  target  di scr im inabi 1 i ty  will  be  possible  until 
theoretical  studies  are  performed  to  parallel  present 
experimental  studies  to  increase  our  understanding  of  these 
factors  discussed  above.  Only  such  studies  will  make 
possible  the  collection  and  analysi-s  of  digital 
multispectral  radiance  data  so  as  to  minimize  errors  of 
omission  and  comission  in  automated  or  semi-automated 
ana  lysis . 


d.  StflSORS 


Sensing  devices  can  consist  of  either  scanners, 
pushbroom  array  devices  or  staring  arrays  129-203,  261-2S6 
The  scanner  has  the  advantage  of  utilizing  a  limited  num¬ 
ber  of  detectors  whose  spectral  response  and  point  spread 
function  may  be  relatively  well  known  (although  this  is 
apparently  not  the  case  for  (e.g.)  the  Landsat  Multi-. 


soectrul  scanner  (MSS).  Neither  the  MSS.  nor  the  TM  had 
documented  inherent  polarization  measurements  performed, 
nor  was  the  spectral  response  o:  each  of  the  channels  or 

each  of  the  bands  on  these  devices  well  documented  in  the 

i .  .  59-62. 

open  literature.  (e.g.  ) 

As  discussed  in  section  5,  tne  point  spread  function  or 
u  sensor  means  that  the  central  portion  of  tne  (  xr  17,  :..i2  a 
higher  sensitivity  than  the  outermost  portions  of  the  IFb'Y. 
For  a  heterogeneous  field  of  view  or  for  a  target  which  does 
not  completely  fill  the  field  of  view,  the  position  of  a 
small  (sub-pixel  sized)  target  with  respect  to  the  peak  of 
tne  point  spread  function  is  therefore  obviously  important. 
This  fact  appears  so  far  to  have  received  relatively  little 
attention  and  it  is  recommended  that  more  work  is  needed 
here.  Pushbroom  scanners  may  consist  of  a  single  array  of 
sensors  which  obtains  information  in  one  spectral  bandpass 
or  may  consist  of  a  two-dimensional  detector  array  placed 
behind  a  wedge  interference  filter  where  the  axis 
perpendicular  to  the  ground  track  obtains  spatially  varying 
information  from  different  pixels  and  the  axis  alongtrack 
obtains  information  in  different  wavebands.  Such  a  device 
clearly  entails  considerable  problems  in  so  far  as  high  data 

rates  involved,  coupled  with  varying  instrument  spectral 
responses  between  detectors,  varying  filter  transmissions 
in  front  of  each  of  the  sensing  elements,  absolute  detec¬ 
tor  calibration  variation  from  element- to-el enent  and 


nheric  :r . tor.  ..’as  allowed  to  varv  from  C.6  no  1.00  in  in 

crements  of  0.01.  For  band  2  (11.0  -  12  am)  the  atmos¬ 
pheric  transmission  was  allowed  to  vary  from  0.60  to  ^.00 
in  incre...ur. ts  of  0.04.  This  resulted  in  a  range  of  ratios 
( t o / 1 ^ ) •  Tnw  coefficients  of  variation  in  both  ana  t0 

was  allowed  to  vary  from  0.01  to  0.05  in  0.01  increments 
which  was  an  estimate  based  upon  the  literature. 

The  discrimination  function  (D.F.)  described  in  equa¬ 
tion  (5.19)  was  plotted  as  a  function  of  the  relative  at¬ 
mospheric  transmissions  and  ~ ^  in  the  two  bandpasses 
used  and  as  a  function  of  the  temperature  of  the  target 
T , .  The  target  temperature  was  allowed  to  vary  from  265° 

L<  to  285°  K  in  one  degree  increments.  The  background  tem¬ 
perature  was  always  275°  K.  The  target  was  assumed  to 
have  a  spectral  invariant  emissivity.  Examples  of  the 
three-dimensional  plots  are  shown  in  Figs.  Cl  to  C3.  The 
vertical  axis  is  DF,  while  the  horizontal  axes  are  the 
ratio  of  the  atmospheric  transmissions  in  the  two 

bands  used  and  the  target  temperature  T^.  Fig.  Cl  shows 
the  case  where  the  correlation  coefficients  between  the 
recorded  radiances  is  0.01,  while  the  coefficient  of  variation 
in  each  of  the  atmospheric  transmission  values  :  .  and  : 


for  bands  1  and  2  is  0.01.  rig.  C2  shows  a  case  which  is 
identical,  'except  that  the  coefficients  of  variation  for 
: ^  and  are  0.05:  clearly  the  increased  variation  in 

atmospheric  transmission  reduces  the  discr iminabil ity 
(i.e.  increases  D.F.)  for  any  given  temperature.  The 
D.F.  is  relative!’.’  insensitive  to  t0/  :  ,  in  either  case. 
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must 

be  less 

(5.18) 


tnan  1.00  for  discrimination  of  target  from  background  using 
the  recorded  radiance  ratios  ^A=  /L^  and  Rg  =  Lg  /Lg^ 
the  95%  confidence  level  (equation  5.15)  Is  given  by: 


at 


D.F. 


(5.19) 


Calculations  of  D.F.  were  made  and  tabulated.  Here,  as  for  the  single 
band  calculations,  samples  of  20  pixels  were  assumed  to  exist  in  the 
target  and  in  the  background  areas.  For  bandl  (3.5-A.ipm)  the  atmos- 


Starting  irom  this  point,  we  may  use  arguments  similar 
to  those  embodied  in  equations  (5.7)  and  (5-8)  to  deduce 
that  the  major  contribution  to  the  variance  in  the  ratio 
between  bands  1  and  2  of  the  radiance  from  target  A  may  be 
attributed  to  atmospheric  transmission  variation. 


That  is, 


(5.15A) 

2p  1™) 

.(  38l/s 

’  S  \  ^ 

W  l  V2 

iiTi ! 
i 

(5.16) 


Ea2 

since  the  variation  in  e.  and  in  e.  is  such  that  — -  is 

A2  A1  £Ai 

always  constant  and  since  BB^  and  BB^j  are  physical  constants 

which  depend  only  on  blackbody  temperature.  We  may  obtain  an 
estimate  of  the  sample  variance  in  the  radiance  ratio  for  target  A: 


modelling  calculations  were  performed  on 


r  ur  the  r 

ratio  of  the  radiance  from  the  target,  compared  to  that  from 
one  background  in  two  different  bands.  That  is,  the  ratio: 


(5.15) 

:  1 


where  L.  is  the  radiant  emittance  from  target  A  in  band  2 
a2 

(11.0-12.0  yra)  and  L.  is  the  radiant  emittance  from 

1 

target  A  in  band  1  (3. 5-4.1  ym).  B12  is  the  blackbody 
radiant  exitance  from  target  A  in  band  2  while  Bj,  is  the 
blackbody  radiant  exitance  from  target  A  in  band  1. 


If  c  \  and  e,  are  the  average  emissivities  of  the 


'A, 


target  A  in  bands  1  and  2,  then,  from  tne  work  of  (e.g.) 
Maxwell"^02,  it  is  reasonable  to  assume  that  there  is  little 
■wavelength  dependence  of  the  emissivities,  so  that; 


1.0  (approximately) 


e  a 


urther  variations  in  e.  will  be  very  strcrgiy  correlated 

Ao 

l 


to  those  in  ,  so  that  the  ratio; 


'A. 


ll 


will  change  very  little  with  variations  in  £^:  the  same  rea- 
ning  mav  be  applied  to  the  emissivitv  of  target  B. 


Three  graphs  are  shown  of  the  sam<_  calculations  at  different 
rotations,  in  order  to  best  present  detail.  Like  calculations 
were  performed  for  the  11.0-12.0  uni  bandpass.  Examples  of 
these  calculations  are  shown  in  Figs.  B1  to  B6.  Figures  B1 
through  B3  show  calculations  of  D.F.  as  a  function  of  (T^, 

(1  -  c . ) )  at  different  rotations  about  the  vertical  (D.F.) 
axis  where  x  =  1.00  and  £g  =  0.90.  Figs.  BA  through  B6  show 
similar  calculations,  with  like  rotations  about  the  D.F.  axis, 
for  i  =  0.60  and  =  0.85.  The  D.F.  peaks  where  target  and 
background  emissivities  are  similar,  especially  at  target  tem¬ 
peratures  close  to  that  of  the  background,  as  one  would  expect. 
However,  discrimination  of  target  from  background  with  even 
turbid,  variable  atmospheric  transmission  and  with  variable 
target  and  background  emissivities  appears  possible,  for  some 
circumstances  even  with  temperature  differences  as  small  as  1 
degree  Kelvin.  The  three-dimensional  plots  show  the  inter¬ 
active  effect  of  target- to-background  emissivity  contrast  and 
temperature  differences,  for  different  atmospheric  transmissions. 
As  expected,  the  higher  atmospheric  turbidity  and  higher 

atmospheric  variability  yielded  a  larger  D.F..  Where  the 
temperatures  and  emissivities  of  target  and  background  were 

similar  discrimination  was  impossible.  The  variation  of  the 
recorded  exitance  would  cause  blurring  of  the  characteristic 
radiant  exitance  vs.  wavelength  curves  for  both  target  and 
background.  The  closer  the  emissivity  values,  the  closer 
(and  less  seperable)  the  curves. 


In  calculation.,  performed,  the  covariance  term  was 
always  small.  The  target  temperature  was  allowed  to  vary  at 
1  uegree  K  increments  from  265K  to  285 K  .  The  backgrouna 
temperature  was  always  taken  to  be  275K  .  The  target 

emissivity  was  allowed  to  vary  at  0.10  increments  from  0.05 
to  0.95*  The  emissivity  of  tne  oac kg round  was  taken  to 

be  0.35  ana  (in  a  separate  run'  0.90.  Atmospheric 
transmission  was  taKen  to  be  0.60  to  1.00  in  0.10 
increments.  The  coefficient  of  variation  (standard 

deviation  divided  by  the  mean)  for  the  emissivities  was 
taken  to  be  0.05,  while  that  of  the  atmospheric  transmission 
was  taken  to  be  0.05  also.  Selected  examples  of  the 
results  of  these  comprehensive  calculations  for  the  bandpass 
3. 5-4.1  tm  are  shown  in  Figs.  Al  to  A6 .  For  each  data  plot, 
the  discrimination  function  (D.F.)  is  shown  as  the  vertical 
axis,  while  the  two  horizontal;  mutually  orthogonal  axes 
are  target  temperature  ( T^) ( temperature  of  background  (T^) 
always  assumed  equal  to  275 K  )  and  (1  -  c^) .  Figs.  Al 
through  A3  show,  for  different  rotations  about  the  vertical 
(D.F.)  axis,  the  value  of  D.F.  plotted  as  a  functi  v,  of 
(T,t  (1  -  For  these  calculations  zu  =  90  and  atm.os- 

ri  A  D 

piier ic  transmission  t  =  1.00.  Figs.  A'-  through  A6  show 

similar  calculations,  with  like  rotations  about  the  D.F. 
axis,  for  a  =  0.70  and  =  0.85.  Each  plot  of  D.F.  as  a 
function  of  (T^,  (1  -  t^) )  is  plotted  for  selected  constant 

background  emissivity  and  atmospheric  transmission  values. 


Now,  as  discussed  in  section  3>  the  half-width  of  the 
95  percent  confidence  interval  {or  of  any  other  pre-selected 
confidence  interval)  on  the  mean  radiance  difference  from 
samples  of  pixels  taken  over  target  and  background  must  be 
greater  than  the  difference  in  the  mean  radiance  levels  (for 
the  bandpass  used). 


i .  e . 


(5.13) 


Thus  the  discrimination  function  (D.F.)  may  be  calculated 
from  equations  (5.11)  -  (5.14). 


Thus,  for  a  single  bandpass: 


from  equations  (5-5)  -  (5-10),  the  sample  variances  in  the 
radiance  recorded  from  target  and  from  background  in  a 
selected  band  are,  respectively: 


while  that  from  the  background  was: 


LBr  eBR  1  t r  .  BB2r  <5-6> 

where : 

£Ar  =  emissivicy  °f  c^e  target  A  in  bandpass  r. 

£gr  =  emissivitv  of  the  background  B  in  bandpass  r. 
x  =  atmospheric  transmission  in  band  r. 

BB^r  =  blackbody  emissivity  for  the  target  (1)  at  temperature 

BB2r  =  blackbody  emissivity  for  the  background  (2)  at  temperature  T2 

It  may  be  shown  (e.g.  ^8  )  that  where  a  dependent 

variable  (v)  is  a  function  of  two  other  variables  (x,y), 
then  the  sample  variance  in  the  dependent  variable  (v)  is 
given  by: 


terms 


(5.7) 

Here  the  terms  of  a  lower  order  of  magnitude  are  generally 
negligible  where  the  coefficient  of  variation  of  each  of  the 
variables  •<  ,  y  is  less  than  about  0.20  and  the  covariance  may 

be  expressed  using  the  well-xnown  formula: 

Sxy  =  p  •  (Sxx)^  '  (Syy)^  (5.8) 

where  p  is  the  correlation  coefficient  between  the  two 
sample  distributions. 
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(14) 


showed  that 


atmospheric  transmission  is  maximum  in  the  regions  3-5-4. 1 
urn  and  10.0-12.0  pm.  For  the  purposes  of  this  study,  we 
have,  for  reasons  of  time,  restricted  our  attention  to 

target  ai scr iminabi  1  i ty  in  each  of  tne  banas  3  -  5  —  4  -  0  -m  and 
11.0-12.0  pm  and  for  the  ratio  of  the  radiance  recorded  in 
these  bands.  We  considered  the  effect  of  random  variations 
in  target  and  background  emissivities  and  in  atmospheric 
transmission  on  the  discriminability  of  target  from 
background.  The  case  of  the  thermal  infrared  region  may  be 
seen  from  the  general  principles  referred  to  in  section  3  to 
be  simpler  than  the  optical-reflective  region:  scattering 
processes  cause  fluctuations  in  target  illumination,  with  a 
consequent  effect  on  target  radiance  fluctuation  only  in  the 
optical-reflective  region. 

We  considered  first  the  case  where  a  target  was  to  be 
discriminated  from  background  in  each  of  the  above-mentioned 


two  banaoasses. 


Here  the  radiance  from  tne  target 


was: 


This  expression  may  be  written  as: 


(5.2) 


where  C-^  =  2tt1iC2 

and  C2  =  h_c 
k 

The  total  (i.e.,  integrated  over  the  electromagnetic 
spectrum)  radiant  exitance  of  a  blackbody  is  given  by: 

C1  A 

=  ~  •  “77"  •  T  (Wm-2)  (5.3) 

rA  15 


Of  course,  the  quantities  used  in  the  above  expressions 
could  just  as  easily  be  expressed  in  c.g.s.  units,  so  that  MA 
and  M-poi’  would  appear  in  c.g.s.  units. 

The  in-band  radiance  is  approximately  expressed  as: 


M 


'A  A 


+  3  ( n  A )  2 


+  6nA  +  6 


e 


(5. A) 


This  expression  was  used  to  calculate  tabular  values  of 
blackbody  exitance  for  different  temperatures  and  for 


In  this  case,  because  of  the  broad  aims  of  the  study, 

it  has  been  decided  to  use  the  assumption  that  tne  target 

and  background  act  as  Lambertian  greybodies  with  almost 

wavelength -independent  emissivities  at  any  given  wavelength. 
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(e.g.  ).  This  is  a  simplistic  assumption.  However,  for 

first-order  calculations,  in  the  light  of  conversations  with 
several  researchers  in  the  field  and  after  a  perusal  of 
available  unclassified  literature,  this  was  considered  to  be 
a  reasonable  starting  point. 

The  Planck  radiation  law  may  be  used  to  calculate  the 
total  radiant  exitance  of  a  blackbody  at  a  selected 
temperature.  This  is  shown  in  equation  (5.1). 


where  M  \  is  the  energy  radiated  per  unit  wavelength  per 

-20-1' 

second,  per  unit  area  of  the  blackbody  (Wm  A  )  (radiant  exitance) 

h  =  Planck' s constant  (J  sec  ^ ) 
c  =  speed  of  light  (m  sec’1) 
k  =  Boltzmann  constant  (J  °K  ^) 

A  =  wavelength 

T  =  absolute  temperature  (°K) 


accumulation  of  scene  elements,  each  of  which  has  its  own 


angular  anisotropy,  considerable  difficulty  may  arise  in 
consistent  target  detection,  discrimination,  tracking  and 
quantification,  depending  upon  the  location  of  the  target 
with  respect  to  the  maximum  of  the  point  spread  function  of 
each  detector  element,  if  the  target  occupies  less  than  1  or 
2  pixels.  Preliminary  studies  have  shown  that  there  are  pos¬ 
sibilities  of  improving  contrast  enhancement  by  dynamic  viewing 
or  "dithering"  (e.g.  )  the  two  dimensional  array.  At  present 

very  few  studies  have  been  reported  in  the  unclassified  liter¬ 
ature.  There  is  a  real  need  to  be  aware  of  the  interaction  of 
scene  heterogeneity,  atmospheric  transmission  fluctuations, 
obscurants  ^  and  target  aspect  with  the  sensing  device  in 

order  to  be  able  to  define  the  optimum  parametric  envelope  for 
target  discrimination,  detection,  tracking  and  quantification. 
Much  more  work  appears  to  be  needed  here. 


5.  MODELLING  STUDIES 


A  considerable  body  of  theory  exists 
of  the  interaction  of  electromagnetic 

dielectric  media.  Generally,  these 

idealized  (e.g.  ^  ’  302-326^ 


for  the  modelling 
radiation  with 
models  have  been 


varying  point  spread  function  f ror.  cetectcr  en.  r.vnt  to 
detector  element.  However,  the  advantages  of  such  a  system, 
are  that  it  does  not  entail  moving  components  such  as  the 
oscillating  mirror  or  rotating  mirror  utilized  in  scanning 
devices  and  should  therefore  have  enhanced  reliability.  It 
is  suggested  that  the  success  of  such  devices  will  rely 
heavily  upon  the  development  of  adequate  onboard  preproces¬ 
sing  facilities  (smart  sensor  devices-e.g. 


A  staring  array  consists  of  a  two-dimensional  matrix  of 
either  quantum  detectors  (such  as  silicon)  or  of  pyro¬ 
electric  detectors  which  record  information  in  the  optical 
reflective,  mid-infrared  or  thermal  infrared  parts  of  tne 
spectrum.  In  the  case  where  charge  coupled  devices  (CCD's) 
or  charge  injection  devices  (ClD's)  are  utilized,  most  of 
the  area  occupied  by  the  array  can  consist  purely  of  sensing 
elements.  Where  the  photodetectors  and  the  circuitry 
involved  in  transferring  the  charge  to  different  registers 
is  located  on  a  single  plane,  then  perhaps  only  50^  of  the 
"real  estate"  of  the  array  element  may  be  occupied  by 
sensors.  The  disadvantages  of  the  two-dimensional  array  are 
the  same  as  mentioned  for  the  multispectral  pushbroom 
scanner.  The  point-spread  function  may  look  something  like 
that  shown  in  Figure  13  for  each  element  of  the  array. 
There  may  be  slight  differences  from  element  to  element,  and 
there  may  be  slight  differences  in  the  absolute  calibration 
(  radiome tr ica 1 ly )  of  each  element.  Thus,  when  a  staring 


.v  ray  v  i  a  w  s  .a 


;onsis  ting  of  a 


he  terogeneous 


Fig.  C3  shows  a  case  where  che  coef f iciencs  of  variation 
for  and  are  0-05  and  where  che  correlation  coefficient 
between  radiance  recorded  in  the  two  bands  is  0.90:  clearly  , 
an  increased  correlation  between  radiance  recorded  in  bands 
1  and  2  (3.5  -  4.1  and  11.0  -  12.0  pm)  increases  the  dis- 
criminability  of  target  from  background  (reduces  D.F.).  In 
general  we  found  that  as  might  be  expected,  the  more  variable 
the  atmosphere,  the  larger  D.F..  However,  the  discrimination 
function  (D.F.)  is  generally  <  1.0  where  the  target  temper¬ 
ature  is  within  1°  K  of  that  of  the  background. 

It  is  significant  that  discrimination  is  independent  of 
target/background  emissivity  contrast  using  the  ratio  method 
and  that  for  target /background  contrasts  of  1°K  (or  less) 
discrimination  appears  more  probable  at  the  95 ”L  confidence 
level  than  for  either  of  the  single  bands,  for  the  levels 
of  atmospheric  turbidity  variance  considered.  Thus,  the 
ratio  method  would  appear  superior  to  the  single  -  channel 
method  and  would  seem  to  be  free  of  the  dependence  on 
emissivitv  variation. 


RECOMMENDATIONS  AND  SUGGESTIONS  FOR  FURTHER 


o  .  COR CLU c  -  U.'i , 
WORK 


Here,  we  present  the  conclusions  drawn  from  this  study, 
plus  general  suggestions  which  one  of  us  (MJD)  has  long 
thought  necessary  steps  in  a  research  program  to  increase  the 
accuracy  of  target  detection,  classification,  quantification 
and  tracking. 


1.  The  discrimination  of  targets  from  background  under 
the  selected  atmospheric  conditions  used  in  the  modelling 
calculations  demonstrates  the  feasibility  of  detection  for 
small  (e.g.  1  degree  Kelvin  =  1°C)  temperature  differences. 

It  appears  that,  for  atmospheric  transmission  and  grounc 
emissivity  random  variations  which  may  be  expected  the 
technique  of  using  ratioed  radiance  values  for  the  bandpasses 
3. 5-4.1  urn  and  11.0  -  12.0  cm  for  reasons  discussed  above 
offers  advantages  over  the  discrimination  possible  using 
radiance  data  recorded  in  either  of  the  single  bands  3.5  - 
•.  .1  ..  m  and  11.0  -  12.0  ..m. 


2.  Real  time  contrast  enhan 
stretching,  geometric  stretching 
enhancement  techniques  with  micro-chip 
be  a  possibility,  providing  that  the 
portion  of  the  instantaneous  field  of 
che  maximum  point  spread  function:  (an 
research).  Clear',  y  if  the  target  is 


cement  using  histogram 
or  other  contrast 
technology  appears  to 
target  occupies  that 
view  (IFOV)  which  has 
area  which  needs 
smaller  than  the  pixel 


(IFOV  at  ground  level)  then  the  contrast  between  the  pixel  containing 

the  target  and  the  pixels  containing  only  background 
material  will  depend  upon  the  position  of  the  point  spread 

function  with  respect  to  the  target.  This  is  an  area  which 
needs  considerable  study.  It  is  suggested  that  little  data 
exists  on  the  point  spread  function  of  any  detecting  device 
and  that  little  is  known  about  improvements  in  contrast 
between  target  containing  pixels  and  background  pixels  using 
dynamic  viewing  or  "dithering"  techniques.  This  area  is  a 
major  gap  in  present  technology  and  must  be  closed, 
otherwise  predictable  discrimination,  identification  and 
quantification  of  targets  may  not  be  possible. 

3.  The  inherent  polarization  of  detection  devices  is  a 
relatively  unexplored  field.  Work  in  .the  optical  reflective 
region  (e.g.  ^7  301^  S^QWS  target  radiance  to  be  polarized: 

studies  performed  in  1974  by  Maxwell  et  al  in  the  emitted 
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infrared  part  of  the  spectrum  show  polarization  to 

exist  in  an  angularly  anisotropic  manner  in  the  thermal 
infrared  region.  Although  considerable  angular  anisotropy 
in  polarization  differences  with  vertical  and  horizontal 
polarizations  have  been  detected,  little  is  apparently  known 
about  the  interaction  of  such  polarized  scene  radiance  with 
inherent  sensor  polarization  effects:  indeed  inherent  polarization 
of  detectors  appears  to  have  received  little  attention. 


Quantitative  measurements  are  necessary  here  in  oraer  to 
ensure  predictable  accuracy  of  targeting  data.  It  is 
essential  that  the  interaction  of  polarized  scene  radiance, 
the  depolarizing  effect  of  the  atmosphere,  the  polarized  sky 
radiance  and  the  inherent  polarization  of  the  sensing 
devices  be  investigated  to  provide  optimum  targeting  data. 
That  is,  it  is  essential  to  minimize  failure  to  detect 
targets, to  minimize  misclassif ication  of  targets  and  to 
utilize  all  potentially  available  information  to  detect, 
quantify  and  track  targets. 

4.  The  possible  combination  of  multiband  approaches  in 
which  polarization  effects  are  utilized  appears  attractive, 
since  polarization  differences  in  scene  radiance  or  in 
target  versus  sky  radiance  are  markedly  greater  in  the 
ultraviolet  part  of  the  spectrum  as  compared  to,  for 
example,  the  emitted  infrared  part  of  the  spectrum.  Such 
considerations  are,  it  is  suggested,  of  vital  importance  in 
horizon  and  nadir-viewing  sensors  based  on  satellite  or  high 
altitude  airborne  platforms.  They  may  also  be  of  use  in 
skytracking  camera  and  in  missile  terminal  guidance  systems. 

5-  The  utilization  of  staring  sensors  to  detect  target 
motion  necessitates  the  knowledge  of  the  point-spread 
function  of  each  element  of  the  array  or  of  any  other  device 
which  is  used  to  image  a  scene.  This  is  because  a  movement 
of  '  less  than  1  pixel  in  sensor  pointing  will  change  the 
radiometric  and  spectral  signature  from  the  scene  because  of 


the  movement  of  the  point-spread  function  with  respect  to 
the  target/background  composite  scene. 


6. 


The  effects  of  scene  obscurants  and  their 


interaction  with  elements  of  two-dimensional  arrays  needs  to 
be  further  studied.  Clearly,  each  element  will  have  a 
slightly  different  absolute  radiometric  sensitivity, 
spectral  response,  inherent  polarization  (especially  in  the 
case  of  ferroelectric  pyroelectric  devices)  and  also  point 
spread  function.  Such  inter-detector  differences  can  lead  to 
errors  in  an  image. 

7.  Cirrus  and  unresolved  cloud,  as  well  as  localized 
obscurants  will  manifest  themselves  as  scene  elements  (i.e. 
components  of  a  heterogeneous  pixel).  The  effects  of  these 
components  upon  radiometric  signatures  needs  considerable 
study  in  so  far  as  it  effects  target  detection,  tracking  and 
quantification.  Modelling  studies  would  be  inexpensive  and 
a  logical  first  step. 

8.  A  model  for  the  two-dimensional  array  in  which  the 
point-spread  functions  are  mapped  in  a  two-dimensional 
diagram  is  necessary  before  it  is  possible  to  determine  the 
optimum  method  of  dynamic  viewing.  It  is  suspected  that 
this  will  depend  upon  the  angular  viewing  regime  and  upon 
the  target- to  -  background  contrast  as  well  as  upon  the 
portion  of  the  pixel  occupied  by  the  target  and  upon  target 
aspect . 


9.  Further  modelling  studies  should  be  performed  and 
should  integrate  terrestrial  vegetation  emittance  models  in 
order  to  predict  parametric  envelopes  for  optimum  data 


acquisition  and  analysis.  For  reasons  discussed  in  section 
3,  it  will  be  necessary  to  model  sensor  output  for  each 
sensor  element  in  a  two-dimensional  array,  pushbroom  array, 
or  scanning  device  so  as  to  take  into  account  the 
radiometric  point  or  line  spread  functions  spectral 
response,  polarization  characteristics,  detector  noise,  and 
the  interaction  of  each  of  these  parameters  with  spectral 
scene  radiance.  The  importance  of  the  difference  in  angular 
anisotropy  in  scene  radiance  from  each  scene  component  will 
depend  upon  the  point-spread  function  of  each  scene  element. 
Models  of  typical  arrays  could  be  used  mathematically  to 
predict  optimum  dynamic  viewing  conditions  to  optimize  target 
to-background  discrimination  under  different  atmospheric  and 
meteorological  conditions,  for  different  geographical  regions 
and  for  different  sensor  types.  Further,  utilizing  the  tech¬ 
niques  of  modelling  studies  could  suggest  possible  new  sensor 
combinations  to  optimize  target  discrimination,  tracking  and 
quantification . 


10.  It  is  strongly  recommended  that  this  work  be 
continued  and  that  modelling  studies  incorporating  some  or 
all  of  the  above  be  pursued  as  a  matter  of  urgency  in  order 


to  better  define  optimum  parametric  envelopes  for  optimal 


data  acquisition  and  real-time  processing.  It  has  been 
shown  that  targets  and  backgrounds  which  are  similar  in 
temperature  may  be  discriminated  under  a  variety  of 
atmospheric  conditions.  The  suite  of  data  used  in  such 
studies  should  be  expanded,  the  models  made  more 
sophisticated  so  as  to  include  spectral  and  point-spread 
function  information  on  the  sensing  devices.  Random  and 
systematic  error  sources  described  above  should  also  be 
further  investigated  in  order  to  determine  their  effects  on 
target  detection,  tracking  and  quantification  accuracy. 


Acknowledgement 


Appreciation  is  expressed  for  support  of  this  work 
under  U.S.  Air  Force  contract  F  30602-81 -C-01 93 •  We 
appreciate  the  assistance  of  D.  Piwinski  with  computing  and 
extensive  work  by  Melanie  Johnson  in  library  research,  word 
processing,  collation  of  data  and  deciphering  the 
handwriting  of  the  principal  investigator. 


60 


FIGURE  CAPTIONS 

Fig.  1.  Overall  picture  of  interactions  between  incident 
and  reflected  radiance,  the  atmosphere  and  the 
earth's  surface. 

Fig.  2.  Interactions  between  solar  radiance  and  the  earth's 
atmosphere . 

Fig.  3.  Nomenclature  for  the  interactions  of  light  energy 
at  the  earth's  surface.  Hemispherica 1 -d i rec tiona 1 
reflectance . 

Fig.  4.  Angular  notation  relating  incident  and  reflected 
radiance.  Bidirectional  reflectance . 

Fig.  5-  Angular  notation  for  reflected  target  radiance. 

Fig.  6.  Wavelength  dependence  and  view  zenith  angle 

dependence  of  bidirectional  reflectance  factor  for 

wheat  at  boot  stage,  for  solar  zenith  and  azimuth 

angles  typical  for  NOAA-6  ephemeris.  The  LARSPEC 

36,70 

data  base  was  used  in  this  calculation 

Fig.  7.  The  relationship  between  incident,  absorbed  and 
reflected  radiance. 

* 

Fig.  8.  The  relationship  between  radiance  from  the  target, 
absorbed  radiance  and  radiance  scattered  from  the 
beam  to  the  sensor. 


Fig.  9-  Wavelength  dependence  of  those  factors  which 
determine  sensor  output. 

Fig.  10.  Showing  that  two  sensors  with  the  same  nominal 
half-power  and  zero-power  bandpass  values  and 
wavelength  limits  can  differ  so  as  to  produce 
different  sensor  outputs  when  viewing  the  same 
target. 

Fig.  11.  Flow  chart  showing  those  factors  controlling 
variations  in  sensor  output  in  the  optical 
reflective  part  of  the  spectrum. 

Fig.  12.  Flow  chart  showing  those  factors  controlling 

variations  in  sensor  output  in  the  thermal  infrared 
part  of  the  spectrum. 

Fig.  13.  The  superposition  of  a  hypothetical  contoured 
point-spread  function  for  a  detector  on  a 
heterogeneous  pixel,  showing  that  the  movement  of 
the  instanteous  field  of  view  by  a  fraction  of  a 
pixel  can  change  the  sensor  output. 

Fig.  14.  The  dependence  of  the  vegetation  index  contrast 

ratio  ( VIN  max MIN  )  on  zenith  angle  and  on 

the  percentage  of  unresolved  cloud  in  the  IFOV  over 

* 

.  the  targets  to  be  distinguished  (percentage  of 

cloud  assumed  to  be  the  same  in  each  case).  Here 
(VIN  =  AVHRR2/AVHRR1  )  and  the  atmosphere  is  clear 
(meteorological  range  >  50  km). 


Fig.  15.  The  dependence  of  the  contrast  ratio  { V  IN  VIN 
on  view  zenitn  angle  and  on  the  percentage  of 
unresolved  cloud  in  tne  IFOV  over  the  targets  to  be 
distinguished  (percentage  of  cloud  assumed  to  be 
the  same  in  each  case).  Here  (VIN  =  AVHRR2/ AVHRR1  ) 
and  the  atmosphere  is  turbid  (meteorological  range  < 
10  km ) . 

Fig.  16.  The  relationship  between  emitted  thermal  target 
radiance,  atmospheric  extinction  and  radiance 
incident  on  the  sensor. 

Figs.  A1-A3.  The  calculated  relationship  between  the  dis- 
criminability  function  (D.F.)  for  detector 
bandpass  3. 5-4.1  cm,  target  temperature  TA 
and  (1  -  e.)  where  c,  *  target  emissivity 
for  background  temperature  Tg  =  275°K,  back¬ 
ground  emissivity  eg  =  0.90  and  atmospheric 
transmission  x  =  1.00.  The  degree  of  varia¬ 
tion  in  emissivity  and  atmospheric  transmis¬ 
sion  considered  is  discussed  in  the  text. 

Three  different  rotations  of  the  same  calcu¬ 
lated  data  are  shown. 

Figs.  A4-A6 .  Similar  to  Figs.  A1-A3,  but  for  x=  0.70  and 
£  g  =  0.85. 

Figs.  B1-B3.  Similar  to  Figs.  A1-A3  but  for  detector  bandpass 
11.0  -  12.0  um:  x  =  1.00  and  eg  =  0.90. 

Figs.  B4-B6.  Similar  to  Figs.  B1-B3  but  for  x  =  0.60  and 
eD  =  0.85. 


Fig.  Cl.  The  discrirainabil ity  function  (D.F.)  is  shown  as  a 


function  of  the  ratio  of  the  atmospheric  transmission 
values  t ^  and  r0  in  two  bandpasses  and  as  a  function 
of  target  temperature  T^.  The  correlation  coefficient 
between  the  radiance  recorded  from  the  target  in  each 
bandpass  is  taken  to  be  0.1  Here  the  coefficients 
of  variation  in  the  emissivities  of  target  r_  .  and 

""  C\ 

of  background  and  in  atmospheric  transmission 
values  t j  >  t 2  are  taken  as  0.01.  The  emissivity 
of  the  target  is  assumed  to  be  spectrally  invariant. 

Fig.  C2.  As  Fig.  Cl  but  with  the  coefficients  of  variation  in 
eA’  eB’  T1  an<^  ~2  =  ®-05. 

Fig.  C3.  As  Fig.  Cl  .  but  with  the  correlation  coefficient 


between  the  radiance  values  recorded  from  the  tar¬ 
get  in  bands  1  and  2  (3.5  -  4.1  ym  and  11.0  -  12.0  pm 
respectively)  taken  to  be  0.90. 
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